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ABSTRACT 


^Incompressible  boundary-layer  transition  has  been  analyzed  using  a 
second-order  closure  turbulence  model.  With  no  transition-specific 
modifications,  the  turbulence  model  predicts  salient  features  of 
incompressible,  zero-pressure-gradient  boundary-layer  transition 
including  sensitivity  to  freestream  turbulence  and  surface  rough¬ 
ness,  transition  width,  and  transitional  velocity  profiles.  With 
transition  modifications  based  on  linear  stability  theory,  the 
model  accurately  predicts  transition  sensitivity  to  surface  heat 
transfer,  pressure  gradient,  and  suction.  With  no  further  modifi¬ 
cations,  transition  predictions  have  been  made  for  several 
hydrodynamic  bodies.  Including  effects  of  surface  heating.-^ 


piTH 
j  crrc  Taa 
I  .j 

■  ’t  t 

I  ^  1  -1  h  t,  ^  ^  * 


il 


ACKNOWLEDGEMENT 


The  authors  gratefully  express  their  appreciation  to  Carl  Gazley  of 
The  RAND  Corporation  for  providing  timely  Information  needed  during 
the  course  of  this  study.  Also,  Joel  Granet  of  NAVSEA  provided 
important  guidance  in  expediting  performance  of  contract  tasks. 


lii 


CONTENTS 


SECTION  PAGE 

ABSTRACT  .  11 

ACKNOWLEDGEMENT  .  Ill 

CONTENTS  .  Iv 

LIST  OP  ILLUSTRATIONS  .  v 

NOTATION  .  vll 

1  INTRODUCTION  .  1 

2  EQUATIONS  OF  MOTION  .  3 

2.1  The  Turbulence/Transltlon  Model  .  3 

2.2  Boundary  Conditions  .  5 

2.3  Transition  Modifications  .  8 

2. H  What  Happened  to  the  Tollmlen-Schllchtlng 

Waves?  .  1^ 

3  BOUNDARY  LAYER  APPLICATIONS  .  17 

3.1  Freestream  Turbulence  .  17 

3.2  Surface  Roughness  .  23 

3.3  Pressure  Gradient  .  23 

3. ^  Suction  .  26 

3.5  Surface  Heat  Transfer  .  29 

^  HYDRODYNAMIC  BODY  APPLICATIONS  .  32 

^.1  Unheated  Bodies  .  32 

4.2  Heated  Bodies  .  36 

5  DISCUSSION  .  40 

APPENDIX:  NEAR-SURFACE  BEHAVIOR  OF  SOLUTIONS 

TO  THE  MODEL  EQUATIONS  .  4l 

REFERENCES  .  45 


LIST  OF  ILLUSTRATIONS 


FIGURE 

1 

2 

3 

il 

5 

6 

7 

8 

9 

10 

11 

12 

13 


Neutral  stability  curves  for  laminar  boundary 
layers  with  pressure  gradient  and  suction  . 

The  stability  function  f(A)  . 

Computed  variation  of  skin  friction  with  Re^  for 
transitional  boundary  layers;  transition  region 
lies  between  points  at  which  minimum  and  maximum 
Cf  occur  . 

Effect  of  freestream  turbulence  Intensity  on 
transition  location  for  a  flat-plate  boundary 
layer  . 

Comparison  of  computed  and  measured  Reynolds 
number  based  on  transition-zone  width  . 

Comparison  of  computed  and  measured  transitional 
velocity  profiles  for  an  Incompressible  flat- 
plate  boundary  layer  . 

Effect  of  surface  roughness  on  flat-plate  boundary 
layer  transition  . 

Comparison  of  computed  and  measured  effects  of 
pressure  gradient  on  boundary-layer  transition  . . . 

Predicted  variation  with  freestream  turbulence 
intensity  of  the  minimum  volume  coefficient 
required  to  delay  transition  on  a  flat-plate 
boundary  layer  with  uniform  suction  . 

CQmin  a  function  of  for  various  freestream 

intensities,  T'  . 

Comparison  of  computed  and  measured  effects  on 
transition  of  neat-transfer  on  low-speed  aero¬ 
dynamic  boundary  layers  . 

Predicted  transition  Reynolds  number  as  a  func¬ 
tion  of  surface  heat  transfer  for  a  low-speed 
aerodynamic  boundary  layer  . 

Body  shapes  and  pressure  distributions  for  the 
four  low-drag  hydrodynamic  bodies  . 


14 


Transition  location  for  the  four  hydrodynamic 
bodies  without  surface  heating  . 


LIST  OF  ILLUSTRATIONS  (concluded) 


FIGURE 


PAGE 


Effect  of  surface  heating  on  Incipient  transition 
Reynolds  number  for  Body  2  . 

Transition  location  for  Body  ^  with  and  without 
surface  heating  . 


NOTATION 


SYMBOL 


e 


e 

max 

r(A),P{A;T;a^} 

h 

H(A) 

k 

L 

e'  0) 

P 

Prj^jPrT 

q 

r 

Rej 


Re,p 

s 

s  s 
min ,  max 

s 

mr 


DEFINITION 
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DEFINITION 

Mean  temperature 

Turbulence  intensity  at  boundary-layer  edge 
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Friction  velocity,  /t  /p 

w  w 

Fluctuating  velocity  in  y  direction 
Body  volume 
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Parameters  in  model  equations 
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1.  INTRODUCTION 


Current  design  techniques  for  low-drag  hydrodynamic  bodies  require 
accurate  prediction  of  transition  from  laminar  to  turbulent  flow. 

For  example,  guided  by  linear  stability  predictions,  judicious  use 
of  various  factors  such  as  surface  heating  and  pressure  gradient 
has  delayed  transition  to  arclength  Reynolds  numbers  in  excess  of 
million  for  small  hydrodynamic  bodies.  Based  on  this  sue  -ss , 
larger  laminar-flow  vehicles  have  been  proposed  for  which  tr-  .sl- 
tlon  Reynolds  numbers  must  be  in  excess  of  200  million.  Agt 
design  procedures  depend  mainly  upon  linear  stability  theory 

Because  of  the  expense  involved  in  building  these  larger  vehicles, 
transition  sensitivity  to  several  additional  factors  must  be  estab¬ 
lished  to  insure  feasibility  of  the  design.  Potentially  detrimen¬ 
tal  effects  of  surface  roughness,  freestream  turbulence,  vibration 
and  acoustic  disturbances  must  be  established.  Because  of  the 
paucity  of  experimental  data  pertaining  to  such  factors,  the 
designer  must  turn  to  theoretically-based  predictive  methods. 

While  linear  stability  theory  is  adequate  for  predicting  effects 
of  primary  design  parameters  such  as  pressure  gradient  and  surface 
heating,  stability  theory  has  no  natural  way  of  simulating  effects 
of  surface  roughness  and  freestream  disturbances.  Alternate 
methods  which  simulate  transition  sensitivity  to  roughness  and 
freestream  disturbances  thus  have  potential  utility  in  the  design 
of  large  laminar-flow  vehicles. 

One  such  method  is  based  on  second-order  closure  turbulence  models. 

As  shown  by  Wilcox  and  Chambers, transition  sensitivity  to  many 
of  the  effects  pertinent  to  hydrodynamic  boundary  layers  can  be 
predicted  with  reasonable  accuracy  using  turbulence  model  equa¬ 
tions.  While  this  method  displays  great  potential  for  engineering 
design,  confidence  in  the  basic  formulation  requires  further  bolster¬ 
ing.  Two  key  points  require  clarification.  Most  significant,  an 
explanation  is  needed  for  the  manner  in  which  Tollmien-Schllchting 
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waves  manifest  themselves  in  the  theory.  The  second  point  is  that 
model  predictions  display  a  strong  sensitivity  to  freestream  dis¬ 
turbances  which  has  been  substantiated  by  limited  experimental 
measurements;  further  substantiation  is  needed. 

The  purpose  of  this  study  has  been  twofold.  First,  an  explanation 
has  been  sought  for  the  way  in  which  the  turbulence-model  approach 
accommodates  Tollmien-Schlichting  waves.  Secori,  turbulence-model 
transition  predictive  accuracy  has  been  assessed  for  the  various 
effects  pertinent  to  transition  on  hydrodynamic  vehicles,  including 
freestream  disturbances.  Section  2  presents  the  model  equations 
upon  which  the  study  is  based.  Appropriate  transition  modifica¬ 
tions  to  the  fully  turbulent  form  of  the  equations  are  devised,  and 
the  model's  relation  to  stability  theory  is  delineated.  Section  3 
summarizes  results  of  classical  incompressible  boundary-layer 
applications  including  comparison  of  predicted  and  measured  effects 
on  boundary-layer  transition  of  freestream  turbilence,  surface 
roughness,  pressure  gradient,  suction  and  surface  heat  transfer. 

In  Section  4,  the  model  is  used  to  predict  transition  on  four  hydro 
dynamic  bodies  with  and  without  surface  heating.  The  concluding 
section  summarizes  results  and  conclusions. 


2.  EQUATIONS  OF  MOTION 


The  turbulence/transition  model  equations  are  summarized  in  this 
section  including  established  values  of  all  closure  coefficients. 
Subsection  2.1  presents  the  model  equations  including  physical 
meanings  of  turbulence  field  properties.  Subsection  2.2  specifies 
surface  boundary  conditions  for  flow  over  rough  surfaces  with  mass 
transfer.  Subsection  2.3  gives  details  of  special  modifications 
needed  to  obtain  accurate  predictions  for  effects  on  transition  of 
pressure  gradient  and  surface  heat  transfer.  In  Subsection  2.4, 
the  manner  in  which  Tollmlen-Schlichting  waves  manifest  themselveb 
in  the  theory  is  explained. 


2.1  THE  TURBULENCE/TRANSITION  MODEL 

Under  the  standard  boundary-layer  approximations,  the  model  equa¬ 
tions  for  two-dimensional  (j=0)  and  axlsymmetrlc  (J=l)  incompres¬ 
sible  flows  are 

Mass  Conservation 

Momentum  Conservation 


Turbulent  Mixing  Enerf 


3  e  .  3  e 


=  [«»l|^|  -  B»u)]pe  +  A|^(^+a»pe)||.] 
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where  s  and  y  are  orthogonal  coor-’lnates  with  s  lying  along  the 
body  (arclength)  and  y  being  normal  to  the  surface;  r  is  the  radial 
distance  from  the  body  axis.  Mean  velocity  components  in  the  s 
and  y  directions  are  denoted  by  u  and  v  while  h  is  the  specific 
enthalpy;  p,  p  and  y  are  mean  density,  pressure,  and  molecular 
viscosity;  t  and  q  are  the  shear  stress  and  normal  heat  flux.  The 
turbulent  mixing  energy,  e,  and  the  turbulent  dissipation  rate,  oj, 
are  needed  to  define  the  eddy  diffuslvity,  e,  which  is  given  by  the 
following  equation: 


e  =  e/u) 

The  shear  stress  and  heat  flux  are 

T  =  (y+Pe)|y 

^  VPrj^  Pr,pj3y 

where  Prj^  and  Pr^p  are  laminar  and  turbulent  Prandtl  numbers, 
quantity  I  is  the  turbulent  length  scale  defined  as 

I  = 


(6) 

(7) 

(8) 


The 


\  / 
e  (1) 


(9) 


The  turbulent  Prandtl  number,  Pr^,  and  the  closure  coefficients 
a,a*,e,3*,a,o*  appearing  in  Equations  (4)  and  (5)  are 


6  =  ^ 


1 

20 


a  = 


_1_ 

2 


3»  = 

^  100 


0.  =  ^ 


Pr  =  — 
T  9 


(10) 


a  = 


a 


•y[l  -  (1  -  A)  exp  (-Re^/2)] 
^1^1  -  (1  -  A)  exp  (-2Re^)J 


(11) 


where  Re^  is  the  turbulent  Reynolds  number  defined  by 

Re^  =  pe^£/y 

Specification  of  the  closure  coefficient  X  Is  deferred  to 
Subsection  2.3- 
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Consistent  with  the  arguments  of  Wilcox  and  Chambers,  ’  the  turbu¬ 
lent  mixing  energy  is  proportional  to  the  kinetic  energy  attending 
the  fluctuation  of  fluid  particles  normal  to  the  plane  of  shear. 
Letting  v'  denote  the  fluctuating  velocity  component  normal  to  the 
shear  plane  (under  the  boundary-layer  approximations,  shear  planes 
are  parallel  to  the  s  direction),  the  turbulent  mixing  energy  is 
given  by 


e 


(12) 


where  <  >  denotes  time  average. 
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The  physical  meaning  of  oj  has  been  discussed  by  Wilcox.  For  incom¬ 
pressible  boundary  layers,  comparison  of  the  limiting  forms  of  the 
model  equations  and  the  exact  Reynolds  stress  equation  very  close 
to  a  solid  boundary  shows  that  u  is  the  rate  at  which  e  is  dissipated 
into  heat,  mean  kinetic  energy  and  other  fluctuation  modes.  For 
incompressible  flows,  Wilcox  deduced  that 


U)  = 


(13) 


where  v =  y/p  is  kinematic  viscosity. 


2.2  BOUNDARY  CONDITIONS 

Boundary  conditions  for  the  model  must  be  specified  at  a  solid 
boundary  (y=0)  and  at  the  boundary-layer  edge  (y=6).  Suitable 

il  q 

boundary  conditions  have  been  devised  by  Wilcox  and  Chambers; 
for  completeness,  the  boundary  conditions  are  summarized  in  this 
subsection . 


For  mean  flow  properties,  conditions  are  the  same  as  for  a  standard 
boundary-layer  computation.  The  velocity  satisfies  the  no-sllp 
condition  at  y=0  while  either  the  surface  temperature,  T^,  or  the 
surface  heat  flux,  q^.  Is  specified.  Hence 


u  =  0  ,  V  = 

T  =  or  9T/3y  =  "  Q^/^^w 


at  y  =  0 


(IM) 


At  the  boundary- layer  edge,  velocity  and  temperature  assume  pre¬ 
scribed  values;  we  thus  write 


u  =  U. 


T  =  T, 


at  y  =  6 


(15) 


Turning  now  to  the  turbulence  parameters  e  and  u,  edge  boundary 
conditions  are  most  conveniently  expressed  in  terms  of  the  free- 
stream  turbulence  intensity,  T' ,  and  the  turbulent  length  scale, 
il=e^/a).  For  transition  applications,  we  generally  use  (see 
Subsection  3*1) 


e  =  |(TV100)2u2 
e  2  e 

a  =  .ooi^ 

e  “ 


y  =  6 


where  a*  =  3/10  Is  the  high  Reynolds  number  (Rerp  >>  1)  limiting  value 
of  a*  and  T’  Is  given  In  per  cent,  i.e., 

T’  s  100 <u'^  +  v'  *  +  w'  ^>Vu^  (17) 

1 3  e  e  e  e 

Finally,  the  surface  boundary  conditions  for  e  and  u  are 

N  \ 

-  e  w  i 

s;  p"  I 

\  (18) 

a*  v„  1 

whicro  k  Is  roughness  height,  u  is  friction  velocity,  and  c^  =  /u^/U^ 
is  skin  friction.  The  quantity  S  is  a  univesal  function  of  surface 
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roughness  and  mass  injection  which  has  been  determined  by  extensive 
study  of  viscous  sublayer  properties;  S  is  given  by 


where 


3  =  /  36  y  ^  {  8  Y 


1  + 


V  /u  >0 

W  T 


(21) 


I  0  ;  v^/u^  <  0 

In  analyzing  the  viscous  sublayer,  the  two  functions  N  and  N  were 
taken  to  be  unity.  For  transition  applications,  and  vary 
with  the  dimensionless  grouping  Cj./(k/0).  Such  a  dependence  has 
been  introduced  to  facilitate  use  of  a  boundary-layer  computation 
for  predicting  transition;  having  =  1  would  require  an  ellip¬ 

tic  Integration  method  to  account  for  local  flow  separation  between 
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roughness  elements  (see  Wilcox  and  Chambers  ’  for  further  details). 

The  functions  N  and  N  have  been  found  by  numerical  experimentation 
e  oj 

their  postulated  dependence  upon  c^/ik/6)  is 


275  c 


u^k/v  <  5 

T  W 


■  -  4 


(22) 


,  k/0  <  275  c 


,  k/0  >  275  c, 


(23) 
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TRANSIT [ON  MODIFICATIONS 


Of  the  various  closure  coefficients,  model-predicted  transit j on  is 
most  sensitive  to  the  coefficient  X  appearing  in  Equations  (10). 

As  argued  by  Wilcox,^  the  value  of  X  can  be  determined  by  demanding 
that  the  model  equations  predict  that  in  a  Blasius  boundary  layer 
turbulent  fluctuations  are  damped  for  Reynolds  numbers  below  the 
linear-stability-theory  minimum-critical  Reynolds  number,  Re^. 
Having  mixing-energy  production,  a*|3u/3y|e,  less  than  mixing 
energy  dissipation,  6*(jje,  Insures  such  damping.  Using  the 
Blasius  velocity  profile  and  the  smooth  wall  w  profile  (l.e., 
u)=20v/6y^  —  see  Appendix),  the  maximum  plate-length  Reynolds  num¬ 
ber,  Re^,  at  which  dissipation  is  greater  than  or  equal  to  produc¬ 
tion  throughout  the  boundary  layer  is 

Re-  =  750/X2  (24) 

The  Reynolds  number  Re-  will  be  equal  to  9*10'*,  the  accepted  value 
of  Re^,  provided 

X  =  1/11  (25) 

With  X  given  by  Equation  (25) j  Wilcox  and  Chambers^  (see  Subsec¬ 
tion  3*1)  have  shown  that  the  model  accurately  simulates  many 
aspects  of  transition  for  an  incompressible  flat-plate  boundary 
layer  (FPBL)  including  transition  sensitivity  to  freestream  turbu¬ 
lence,  transition  width,  and  transitional  velocity  profiles.  How- 

•5 

ever,  previous  experience  with  heated  hydrodynamic  boundary  layers 
shows  that  in  order  to  accurately  predict  effects  of  pressure  grad¬ 
ient  and  surface  heat  transfer.  Equation  (25)  must  be  modified. 

The  remainder  of  this  subsection  is  devoted  to  explaining  the  need 
for  and  specification  of  further  modifications  to  the  closure 
coefficient  X. 

As  noted  above,  the  value  of  X  has  been  fixed  by  demanding  that 
the  linear-stability  minimum-critical  Reynolds  number,  Re^,  for 
the  Blasius  boundary  layer  match  the  corresponding  model-equation 


neutral-stability  Reynolds  number,  Re^.  Demanding  that  Re  =  Re^ 
for  the  Blaslus  boundary  layer  yields  the  value  of  X  given  in 
Equation  (25).  The  model  equations  reasonably  can  be  expected  to 
apply  to  transitional  flows  which  are  insensitive  to  spectral 
effects.  That  is,  the  various  constants  in  the  model  equations 
are  essentially  correlation  coefficients  which  have  been  inte¬ 
grated  over  the  turbulent  spectrum.  Hence,  if  the  stability  dia¬ 
gram  shows  that  a  wide  range  of  wave  numbers,  a,  undergo  amplifi¬ 
cation,  the  spectrum  will  more  closely  resemble  a  fully-turbulent 
spectrum  than  if  only  a  small  range  of  wave  numbers  are  unstable. 

For  example,  the  stability  diagram  for  a  boundary  layer  subjected 
to  a  pressure  gradient  is  shown  in  Figure  1.  For  adverse  pressure 
gradient,  a  finite  range  of  wave  numbers  are  unstable  at  all 
Reynolds  numbers  in  excess  of  (note  that  6*  is  displacement 

thickness).  On  the  basis  of  the  discussion  above,  the  model  would 
be  expected  to  accurately  predict  the  destabilizing  effect  of 
adverse  pressure  gradient.  In  contrast,  the  stability  diagram 
becomes  thinner  with  increasing  favorable  pressure  gradient  so 
that  spectral  effects  become  increasingly  Important,  particularly 
for  small  T’  which  yields  transition  at  large  values  of  Re^*;  the 
model  hence  would  be  expected  to  fare  poorly  for  transitional 
boundary  layers  with  favorable  gradients  (and  small  freestream 
disturbances ) . 

As  shown  by  Wilcox  and  Chambers,  the  original  version  of  the  model 
behaves  Just  as  the  above  discussion  indicates.  Excellent  agree¬ 
ment  between  theory  and  experiment  is  obtained  for  adverse  grad¬ 
ients  while,  for  low  freestream  turbulence  intensities,  the  model 
falls  to  predict  the  strong  stabilizing  effect  of  favorable  grad- 
lent.  To  remove  this  deficiency,  Wilcox  and  Chambers  introduced 
an  empirical  modification  to  Equation  (25).  While  reasonably  good 
agreement  with  measurements  resulted,  the  modification  lacked  rigor. 

'  ri  the  rmei'ent  study,  a  better  approach  to  modifying  X  has  been 
found.  That  is,  the  requirement 
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has  been  extended  to  Include  favorable  pressure  gradients.  Hence, 
for  small  freestream  turbulence  intensity  we  expect  to  have 

X  =  ^  f(A)  as  T’  0  (27 

where  A  is  the  modified  Pohlhausen  pressure  gradient  parameter 


A 


(28 


The  function  f(A)  must  be  determined  by  equating  Re  and  Re^. 

Figure  2  presents  results  based  on  the  Phlhausen  profiles  ,  a  good 
fit  to  the  data  indicate  the  variation  of  X  with  A  is  hence 

^  ^  FF  [-^OAH(A)]  (29 


where  H(A)  is  the  Heaviside  stepfunctlon.  The  limiting  value  of 
f(A)  as  A 0  has  been  obtained  by  a  similar  analysis  of  the  asymp¬ 
totic  laminar  profile  for  a  uniformly-sucked  PPBL. 


Turning  to  heat  transfer,  an  additional  modification  to  X  is  needed 

•3 

Following  Wilcox  and  Chambers,  for  incompressible  aerodynamic 
boundary  layers  with  heat  transfer,  the  neutral  stability  Reynolds 
number  based  on  wall  conditions  is  given  by 


Re 


730 

X" 


(30 


According  to  linear  stability  theory,  the  corresponding  minimum- 

critical  Reynolds  number  Re^  varies  as  follows: 

w 


Rec  •V  (T„/T^)-’ 
w 


(31 


Hence,  since  y -x.  T®'’  for  air,  we  have 


1 - r 


O  INFERRED  FROM  LINEAR- 
STABILITY  ANALYSIS: 
FOHLHAUSEN  PROFILES 

•  INFERRED  FROM  LINEAR- 
STABILITY  ANALYSIS: 


f(A 


Figure  2.  The  stability  function  f(A) 


(32) 


Re  'V/  ( p  /y  ) 
w 

We  thus  postulate  that  for  flows  with  heat  transfer,  X  be  written 
as  follows: 


X  =  P{A  ;  T'  } 

liypg/  max 


(33) 


Finally,  note  that  for  high  freestream  turbulence  Intensities, 
transition  Is  unlikely  to  be  sensitive  to  spectral  effects  regard¬ 
less  of  the  stability  diagram.  That  is,  typical  high  intensity 
freestream  turbulence  will  have  fluctuations  at  all  frequencies 
(wave  numbers).  Therefore,  the  modification  proposed  in  Equa¬ 
tion  (33)  is  strictly  valid  only  as  T*  ^0.  Hence,  to  complete  the 
formulation,  we  introduce  an  exponential  dependence  upon  the 

maximum  disturbance  in  the  boundary  layer,  so  that  the  completed 
formulation  is  as  follows. 


TRANSITION  MODIFICATION 


1  /^w\^ 


F{A  ;  T’  } 
’  max 


r(A)  =  +  fa®”''’  t-'tOAH(A)] 


where 


9 ^  / 3^u\ 

■pw^eVayVw 


and  T'  =  100 
max 


3  ®max^^e 


2.4  WHAT  HAPPENED  TO  THE  TOLLMIEN-SCHLICHTING  WAVES? 

Before  proceeding  to  applications.  It  Is  Instructive  to  discuss 
the  key  aspect  of  the  turbulence-model  transition-prediction 
method  requiring  clarification,  namely,  the  manner  In  which 
Tollmlen-Schllchtlng  waves  manifest  themselves  In  the  theory. 
Because  we  are  solving  the  conventional  long-time  (Reynolds)  aver¬ 
aged  equations  of  motion,  no  wave-llke  motion  can  be  discerned. 
Hence,  the  whole  stability  mechanism  can  only  be  represented 
Implicitly  with  the  model  equations.  The  arguments  presented  In 
Subsection  2.3j  particularly  the  point  that  the  various  closure 
coefficients  are  correlation  coefficients  which  have  been  inte¬ 
grated  over  the  wave-number  spectrum.  Illustrate  the  manner  in 
which  Tollmlen-Schllchtlng  wave-stability  information  Is  Implied. 
Clearly  the  wave-number  spectrum  of  a  given  correlation  coefficient 
associated  with  Initiation  and  early  amplification  of  Tollmlen- 
Schllchtlng  waves  is  likely  to  be  quite  different  from  the  corre¬ 
sponding  spectrum  In  fully  turbulent  flow.  Thus,  the  values  of 
the  closure  coefficients  a  and  a*  at  the  beginning  of  transition 
differ  markedly  from  the  values  appropriate  for  fully  turbulent 
flows . 


Analysis  of  the  viscous  sublayer  indicates  that  a  and  a*  are 
smaller  at  low  Reynolds  numbers  than  at  higher  Reynolds  numbers, 
the  ratio  being  X,  i.e.. 


a»(Re^  =  0) 
a»(Re^  oo) 


(35) 


and  similarly  for  a.  While  the  value  quoted  in  Equation  (25)  Is 
sufficient  to  yield  accurate  sublayer  structure  and  transition 
predictions  for  the  PPBL,  we  have  found  that  a  constant  value  for 
A  Is  insufficient  to  provide  an  accurate  model  for  more  general 
transition  applications.  It  Is  through  the  precise  value  of  X 
that  the  Tollmlen-Schllchtlng  waves  implicitly  appear  in  the 
model.  Because  the  only  transition-specific  modifications  to  the 
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model  equations  are  for  the  coefficient  X,  It  is  hence  unsurprising 
that  ultimately  we  have  chosen  to  rely  upon  linear  stability  theory 
to  set  its  value. 

Given  this  insight,  the  whole  concept  of  using  turbulence-model 
equations  to  describe  transition  can  be  cast  in  a  different  light. 

On  the  one  hand,  time-averaging  conceals  many  Important  phys¬ 
ical  aspects  of  transition  mechanisms,  particularly  during  the 
early  linear-amplification  phase.  Vie  are  thus  obligated  to  put 
some  of  the  physics  back  into  the  equations  which  was  lost  through 
the  time-averaging  process;  ergo,  the  modification  to  X.  On  the 
other  hand,  assuming  the  latter  phases  of  transition  to  be  very 
rapid,  the  time-averaging  process  becomes  a  more  plausible  concept 
as  the  flow  more  nearly  resembles  a  turbulent  flow.  Interestingly, 
the  modification  to  X  has  little  effect  on  the  latter  stages  of 
transition.  The  turbulence-model  transition-prediction  approach 
thus  has  its  strongest  foundation  in  the  latter  stages  of  transi¬ 
tion,  precisely  the  regime  where  conventional  linear-stability 
methods  are  not  well  founded. 

Q 

To  explain  this  last  point,  note  that  the  classical  Smlth-Gamberonl^ 
e’  method  Ignores  nonlinear  effects  in  a  region  where  such  effects 
certainly  have  a  strong  effect  on  the  transition  process.  Thus, 
turbulence-model  equations  might  most  properly  be  viewed  as  a 
plausible  alternative  to  the  e®  method.  That  is,  a  linear  sta¬ 
bility  computation  could  be  performed  up  to,  and  perhaps  a  bit 
beyond,  the  minimum  critical  Reynolds  number.  Results  of  the  sta¬ 
bility  computation  would  define  X.  Then  rather,  than  continuing 
to  solve  eigenvalue  problems  to  determine  amplification  factors  up 
to  the  e®  amplification  point,  the  model  equations  could  be  used 
to  predict  transition  location.  What  we  are  currently  doing  is 
using  a  correlation  of  linear  stability  minimum  critical  Reynolds 
numbers  to  fix  the  value  of  X. 
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All  of  the  comments  above  pertain  to  transition  triggered  by 

disturbances  of  sufficiently  low  amplitude  for  linear  stability 

theory  to  be  relevant.  However,  if  large  amplitude  freestream  or 

roughness-induced  disturbances  are  present,  the  bypass  phenomenon^^ 

may  occur,  in  which  case  Tollmien-Schlichtlng  waves  are  irrelevant. 

Because  the  turbulence/transition  model  accurately  simulates 
4 

bypass,  most  notably  for  roughness-induced  transition,  the  cur¬ 
rent  formulation  enjoys  an  interesting  advantage  over  linear- 
stability  methods  for  flows  with  large  amplitude  disturbances. 


3.  BOUNDARY  LAYER  APPLICATIONS 


To  test  the  model,  we  first  consider  carefully-documented  experl- 
rne'nts  for  conventional  Incompressible  boundary  layers.  This 
section  presents  results  of  five  applications,  all  of  which  are 
pertinent  to  low-drag  hydrodynamic  vehicles.  First,  we  apply 
the  model  to  an  Incompressible  FPBL  Including  comparisons  of  com¬ 
puted  and  measured  transitional  velocity  profiles,  transition 
width  and  transition  sensitivity  to  freestream  turbulence  Intensity 
and  scale.  Then,  we  simulate  the  effects  of  surface  roughness, 
pressure  gradient  and  suction.  In  the  concluding  subsection,  we 
use  the  model  to  analyze  effects  of  surface  heat  transfer  on  aero¬ 
dynamic  boundary  layers. 

3.1  FREESTREAM  turbulence"^ 

As  the  first  step  In  testing  the  model,  we  analyze  various  aspects 
of  model-predicted  Incompressible  FPBL  transition.  As  shown  in 
Figure  3j  starting  from  laminar  flow  at  the  plate  leading  edge, 
the  model  predicts  that  the  skin  friction  initially  matches  the 
Blasius  value.  Then,  depending  upon  the  freestream  turbulence 
Intensity,  T' ,  skin  friction  rapidly  Increases  at  a  critical 
Reynolds  number,  Re^,  ,  and  asymptotically  approaches  the  equlllb- 
rlum  turbulent  value.  Predicted  variation  of  c^,  closely  resembles 
that  observed  when  a  boundary  layer  undergoes  transition  to 
turbulence. 

For  example.  Figure  4  shows  that,  consistent  with  measurements,^^ 
momentum-thickness  Reynolds  number  at  transition,  varies 

almost  linearly  with  T’  for  low-intensity  freestream  turbulence 
(l.o.,  T'  less,  than  1%)  .  (The  criterion  used  to  define  transition 


Some  of  the  results  presented  In  this  subsection  were  obtained 
In  Contract  F4^620-7^-C-00^8 . 
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TURBULENT: 

KARMAN-SCHOENHERR 


Figure  3-  Computed  variation  of  skin  friction  with  Re^  for  transitional  boundary 
layers;  transition  region  lies  between  points  at  which  minimum  and 
maximum  Cf  occur. 


is  the  point  at  which  achieves  a  minimum.)  Note  that  turbulence 
scale  has  an  effect  on  Re©^ ,  particularly  for  T’  >1^.  Computations 
have  been  performed  with  various  ratios  of  /a^  to  6  ;  having 
Z  //a^ 6  given  by 

Z 

=  .004  ;  T’  <  1  (36) 


most  nearly  matches  the  low-intensity  values  of  Re, 


Somewhat 


larger  values  of  are  needed  to  match  the  high-intensity 

data;  for  values  of  T'  in  excess  of  1%,  excellent  agreement  between 

computed  and  measured  Re0  is  obtained  with  given  by 

t  e  “ 


/oiy  6 


=  .004T'  ;  T’  >  1 


Intuitively,  we  expect  that  Z^  should  Increase  with  T'  for  hlgh- 
intenslty  turbulence  since,  in  the  case  of  fully  turbulent  boundary 
layers,  values  of  generally  are  an  order-of-magnltude 

greater  than  that  given  by  Equation  (36).  The  variation  of  Z^  with 
T'  given  by  Equation  (37)  is  quite  plausible  since  the  peak  local 
intensity  in  a  fully  turbulent  incompressible  FPBL  is  of  the  order 
of  10  to  \2%  which,  from  Equation  (37)  indicates  Z  //^6  is  of  the 
order  .04  to  .05.  For  turbulent  boundary  layers,  is 

typically  .09- 


In  Figure  5,  predicted  width  of  the  transition  region  is  compared 
12 

with  measured  width.  Transition  width,  Ax^,  is  defined  as  the 
distance  between  minimum  and  maximum  skin-friction  points.  Using 
this  definition  for  Ax^,  Reynolds  number  based  on  Ax^  was  computed 
for  R^x^  ranging  from  5.0*10‘'  to  4.4*10®.  As  shown,  the  computed 
curve  falls  within  experimental  data  scatter. 

Figure  6  exhibits  computed  and  measured^^  velocity  profiles  through 
transition  for  T'  =  .03%.  Comparison  of  velocity  profiles  at 
x=  5.75  ft  Indicates  the  computed  boundary  layer  goes  turbulent  a 
bit  faster  than  measured.  Proceeding  downstream,  however,  computed 
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Comparison  of  computed  and  measured 
Reynolds  number  based  on  transition 
zone  width. 


and  measured  velocity  profiles  show  decreasing  differences  until, 
beyond  x=6.'75  ft,  differences  are  less  than  5%- 

An  Interesting  point  about  these  results  Is  that  because  (9^u/9y^) 

w 

vanishes  for  Incompressible  FPBL  flow,  the  transition  modification 
defined  In  Equation  (34)  has  no  effect.  Hence,  the  turbulence  model 
as  stated  In  Subsection  2.1  applies,  with  X=l/ll,  to  the  Incompres¬ 
sible  FPBL. 

3.2  SURFACE  ROUGHNESS 

Figure  7  shows  computed  effects  of  surface  roughness  on  FPBL  transl- 

l4 

tlon;  experimental  data  of  Felndt  are  shown  for  comparison. 
Freestream  Intensity  has  only  a  slight  effect  on  transition  Reynolds 
number  based  on  plate  length,  R^x^j  f'or  roughness-height  Reynolds 
numbers  in  excess  of  300.  Thus,  consistent  with  qualitative  obser¬ 
vations,  the  model  predicts  existence  of  a  roughness  dominated 
regime,  although  substantiating  data  are  unavailable  for  determin¬ 
ing  the  minimum  value  of  Rej^  at  which  transition  becomes  roughness 
dominated.  Also,  again  consistent  with  measurements  and  qualita¬ 
tive  observations,  roughness  has  virtually  no  effect  on  FPBL 
transition  for  Rsj^  ^  120. 

3.3  PRESSURE  GRADIENT 

Figure  8  compares  predicted  and  measured  effects  of  pressure  grad- 
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lent  on  FPBL  transition.  As  shown  In  the  figure,  the  Crabtree 

data  are  closely  simulated  by  the  computations  with  T'  =  .01%  and 

.03%.  Since  the  Crabtree  data  are  for  flight  tests  (closed  circles) 

and  quiet  wind  tunnels  (open  circles),  the  value  of  T'  for  the 

experiments  would  be  expected  to  fall  In  the  range  .01%-. 03%.  The 

agreement  between  theory  and  experiment  hence  Is  very  good  for  low 

l4 

Intensities.  Additionally,  the  Felndt  hlgh-lntensity  data 
(T'  =  1.25%)  are  closely  matched  by  the  computed  curve  with  T'  =  1.25% 
because  Felndt  found  that  transition  occurs  ahead  of  the  minimum 
critical  Reynolds  number  for  favorable  gradient,  this  flow  exhibits 
an  example  of  the  bypass  phenomenon. 


3.^ 


SUCTION 


Figure  9  shows  computed  minimum  volume  coefficient,  required 

to  prevent  transition  of  a  FPBL  with  uniform  suction.  Volume 
coefficient  Is  defined  as 


C 


Q 


P  V 
^w  w 

P  u 
e  e 


(38) 


where  v  Is  suction  velocity.  As  shown,  for  values  of  T'  above 
about  0.3%,  .0012,  a  value  typical  ’  of  experiments  per 

formed  In  noisy  wind  tunnels.  As  expected,  our  numerical  computa¬ 
tions  confirm  that  in  the  limit  T'  ->-0  the  model  reproduces  the 
linear-stability  predicted  ^  of  1.4*10”®. 


In  computing  we  have  demanded  that  transition  be  delayed 

Indefinitely,  i.e.,  that 


CO 


(39) 


It  Is  Instructive  to  compute  a  family  of  C 
stringent  condition 


Qmln 


curves  for  the  less 


Re 


>  R. 


min 


(40) 


where  is  a  specified  minimum  Reynolds  number  below  which 

transition  Is  to  be  prevented.  Such  a  family  of  curves  was  con¬ 
structed;  results  are  shown  In  Figure  9*  Consistent  with  the 
1 8 

Lang  data,  provided  T'  is  not  too  large,  of  the  order  of 

.0001-. 0003  Is  sufficient  to  stabilize  a  FPBL  up  to  Reynolds  num¬ 
bers  of  a  few  million.  However,  for  a  fixed  turbulence  level, 

achieving  values  of  Re„.  beyond  10  million  requires  more  suction. 

■''•t 

To  see  this  quantitatively,  is  replotted  In  Figure  10  with 

T’  held  constant  on  each  curve.  As  an  example,  when  T*  =  0.105S, 
a  value  of  =  .0001  is  sufficient  to  delay  transition  to 

Rex^'^6‘10®  while  must  increase  by  a  factor  of  almost  five 

in  order  to  have  Rex4.  =  200*10® . 


Figure  9.  Predicted  variation  with  freestream 
turbulence  Intensity  of  the  minimum 
volume  coefficient  required  to  delay 
transition  on  a  flat-plate  boundary 
layer  with  uniform  suction. 


3.5  SURFACE  HEAT  TRANSFER 

As  our  final  boundary-layer  application  we  now  turn  to  effects  of 

surface  heat  transfer  on  aerodynaunlc  boundary-layer  transition. 

Figure  11  compares  computed  and  measured  transition  Reynolds  number 

for  low-speed  aerodynamic  boundary  layers.  The  data  of  Zyslna- 
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Molozhen  and  Kuznetsova  were  taken  in  relatively  noisy  environ¬ 
ments  so  that  they  correspond  to  high  intensity  freestream  condi¬ 
tions.  The  fact  that  the  curve  computed  with  T'  =  1.25^  is  close  to 
the  data  is  hence  very  encouraging.  For  low  intensity  computations 
(T'  =0.03^),  the  predicted  stabilizing  effect  of  cooling  is  much 
stronger  than  for  the  T'  =  1.25/?  computations.  As  expected,  the 
variation  of  (Rey, )  is  approaching  the  linear  stability  predictions 
ofMackF 

Because  v  increases  with  T  .  a  subtle  feature  of  surface  heating 
w  w 

effects  on  transition  is  masked  by  displaying 

a  function  of  Specifically,  as  the  heating  rate  decreases, 

x^  initially  increases,  achieves  a  maximum,  and  eventually  decreases 
Figure  12  shows  (Re^^^)  U^x^/v^  as  a  function  of  T^/T^  for  the 

T'=1.25%  calculation.  As  shown,  for  T^/Tg<0.5,  additional  cool¬ 
ing  destabilizes  the  boundary  layer.  This  trend  is  consistent  with 
linear  stability  ■predictions. 


HYDRODYNAMIC  BODY  APPLICATIONS 


Having  tested  the  model  for  well-documented  incompressible  flows  on 
simple  flat-plate  geometries,  we  now  turn  to  practical  hydrodynamic 
designs.  Computations  presented  in  this  section  have  been  performed 
with  virtually  no  advance  knowledge  of  either  linear-stability  pre¬ 
dictions  or  ‘.experimental  data.  The  section  consists  of  two  parts. 
First,  we  simulate  transition  on  four  unheated  hydrodynamic  bodies; 
transition  location  is  presented  for  each  body  over  a  range  of 
freestream  velocities.  Then,  for  two  of  the  bodies,  the  stabilizing 
effect  of  surface  heating  on  transition  is  predicted. 


^4.1  UNHEATED  BODIES 


Four  typical,  axlsymmetrlc  low-drag  bodies  have  been  considered. 
Figure  13  depicts  the  body  shapes  and  Douglas-Neumann  computed 
pressure  distributions.  In  the  figure,  L  denotes  body  length  mea- 
sur-ed  along  the  symmetry  axis,  is  the  arclength  at  which  the 

body  radius  is  a  maximum,  and  r^^^  Is  the  maximum  body  radius.  Also, 
Cp  is  the  pressure  coefficient  defined  by 


C 


P 


P-Poo 


(41) 


Computations  have  been  performed  with  molecular  viscosity  and 
laminar  Prandtl  number  given  by  values  appropriate  for  water  near 
room  temperature,-"  viz, 

y  =  7.943‘10"®(T/600)"®  Ib^sec/ft^  (42) 

Prj^  =  2.5l(T/600)"®  (43) 

with  T  in  degrees  Rankine.  Freestream  turbulence  intensity  is 
assumed  to  be 

T’  =  .025S?  (44) 
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CURVE  BODY 


and  SLq  Is  given  by  Equation  (36).  Some  of  the  computations  have 
been  repeated  with  T*  =  .01^  and  =  .09;  transition  Reynolds 

numbers  changed  by  less  than  12??.  Roughness  heights  of  zero, 

16  pin  and  pin  have  been  used;  variation  of  k  over  this  range 
produces  a  negligible  effect  on  transition. 

Inspection  of  Figure  13  shows  that  Bodies  1,  2  and  3  have  progres¬ 
sively  strong  favorable  pressure  gradients  ahead  of  the  maximum 
radius  (minimum  pressure).  Body  A  has  a  mild  favorable  gradient 
up  to  s/s^j,=  0.^1,  followed  by  a  slight  adverse  gradient  from 
s/s^^=0.4l  to  s/s  =0.73j  and  subsequently  a  favorable  pressure 
gradlent  up  to  the  pressure  minimum.  Intuitively,  we  expect  to 
find  Increasing  transition  Reynolds  number  as  we  proceed  from 
Body  1  to  Body  3-  For  Body  4,  we  expect  transition  will  occur 
either  upstream  of  or  very  near  the  beginning  of  the  adverse  pres¬ 
sure  gradient  on  the  forebody. 

Model  predictions  for  the  four  bodies  confirm  our  Intuitive  notions. 
Figure  14  presents  transition  arclength,  s^,  as  a  function  of  volume 
Reynolds  number,  Reyi/s,  defined  by 

Reyi/3  5  pU„V^'Vp„  (45) 

For  Body  1,  transition  occurs  ahead  of  laminar  separation  for  Reyi/s 
In  excess  of  about  2  million.  Because  of  Body  2's  stronger  favor¬ 
able  pressure  gradient,  transition  does  not  move  ahead  of  laminar 
separation  until  ROyi/s  exceeds  9.8  million.  For  the  even  stronger 
favorable  gradient  of  Body  3,  although  transition  appears  just 
upstream  of  laminar  separation  at  Reyi/3  =  6.5  million,  transition 
occurs  downstream  of  minimum  pressure  at  the  highest  Reynolds  num¬ 
ber  considered,  namely  15  million.  Finally,  for  Body  4,  the 
adverse  pressure  gradient  on  the  forebody  triggers  transition  at 
the  lowest  Reynolds  numbers  considered;  for  the  higher  Reynolds 
numbers  transition  occurs  upstream  of  the  adverse  gradient. 

Table  1  summarizes  results  of  the  computations. 
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Transition  location  for  the  four  hydro 
dynamic  bodies  without  surface  heating 


Table  1.  Summary  of  transition  predictions 
for  the  unheated  bodies. 


Body 

Number 

10'"  Re^i/3 

s . /s 
t  mr 

10"®  Re 

^t 

L 

1.12 

2.79 

5.57 

11.17 

1.108* 

0.841 

0.346 

0.178 

4.68* 

8.43 

6.67 

6.78 

2 

5.19 

1.160* 

12.23* 

10.36 

1.070 

22.98 

11.97 

0.847 

20.28 

13.97 

0.661 

17.63 

15.17 

0.568 

16.05 

17.17 

0.471 

14.63 

1 

19.95 

0.378 

13.22 

3 

1.88 

1.253* 

4.14* 

5.63 

1.253* 

12.50* 

9.34 

1.140 

19.13 

13.14 

1.059 

25.39 

15.02 

1.049 

28.75 

3.98 

0.535 

6.55 

7.96 

0.389 

9.61 

11.93 

0.315 

11.60 

15.91 

0.264 

12.77 

*Lamlnar  separation. 


4.2 


HEATED  BODIES 


For  Bodies  2  and  4,  computations  have  been  repeated  with  various 

amounts  of  surface  heating.  Two  heating  distributions  have  been 

considered  for  Body  2.  For  both  distributions  there  Is  no  heating 

up  to  s  =  s  .  ,  uniform  heating  from  s=s.  tos  =  s  ,  and  no 
min  min  max 

mln'^^mr 

=  .17  for  the  second 
distribution.  The  heating  distributions  for  Body  4  are  those  pre¬ 
dicted  by  linear-stability  theory  as  being  the  minimum  heating 
required  to  prevent  transition  from  occurring  ahead  of  laminar 
separation. 


heating  beyond  s  =  s^^^^  (see  insert  In  Figure  15);  s^^^/s^^=  .46 

for  the  first  distribution  while  s  ^„/s 

min  mr 
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l''i)’;ur’c  I'j  and  Table  2  present  results  of  the  computations  for 

Body  2.  Incipient  transition  Reynolds  number,  (Re.-i/s)  ,  (i.e 

v  Inclp 

the  value  at  which  transition  and  laminar  separation  points  coin¬ 
cide)  has  been  predicted  as  a  function  of  heating  rate.  As  shown 

when  s  .  /s  =.46  a  modest  increase  in  (Re.^ws)  is  observed, 

min  mr  V-^'^'lncip 

Moving  closer  to  the  nose  causes  a  more  dramatic  increase  in 

(Re„i/3)  .  This  trend  is  in  qualitative  agreement  with 

V  incip 
measurements . 


Table  2.  Summary  of  transition  predictions  for 
Body  2  with  surface  heating. 


0 

.75 

1.50 


0 

.75 

1.50 


10"®(Re„i/3) 


^  incip 


9.8 

12.0 

14.0 


17.6 


10  ®  Re. 


22.0 

32.0 

40.3 


Figure  16  and  Table  3  summarize  results  for  Body  4  with  surface 
heating.  At  the  two  lowest  Reynolds  numbers,  transition  has  been 
prevented  upstream  of  laminar  separation.  For  Rey  1/3  =  11 . 93  >  transi 
tlon  is  predicted  between  the  minimum  pressure  and  laminar  separa¬ 
tion  points  while,  at  the  highest  Reynolds  number  considered, 
transition  occurs  just  upstream  of  minimum  pressure.  Again, 
predictions  are  in  qualitative  agreement  with  measurements. 


Table  3.  Summary  of  transition  predictions  for 
body  4  with  surface  heating. 


10  *  Rey  1/3 

s./s 
t  mr 

10“®  Re 

3.98 

1.276« 

15.91* 

7.96 

1.276» 

32.10* 

11.93 

1.196 

45.54 

15.91 

0.951 

46.51 

Laminar  separation 


5.  DISCUSSION 


Results  presented  In  Sections  3  and  4  lend  further  confidence  to 
the  turbulence-model  transition-prediction  method.  The  model 
accurately  simulates  transition  sensitivity  to  many  of  the  effects 
pertinent  to  low-drag  hydrodyncimic  bodies.  While  detailed  compari¬ 
sons  between  predicted  and  measured  transition  points  for  the 
hydrodynamic  todies  analyzed  in  Section  4  have  not  been  presented, 
early  Indications  are  that  an  acceptable  level  of  accuracy  has 
been  obtained. 

The  explanation  for  the  manner  in  which  Tollmlen-Schllchting  waves 
manifest  themselves  in  the  theory  (Subsection  2.4)  puts  the  model 
in  much  clearer  perspective  than  previously  realized.  Most  notably, 
by  using  linear  stability  theory  to  determine  the  closure  coeffi¬ 
cient  X  for  a  given  flow  the  model  could  be  the  basis  of  a  plausible 
alternative  to  the  o*  procedure.  Predictions  remain  relatively 
sensitive  to  freestream  turbulence  intensity,  a  sensitivity  which 
probably  will  remain  uncertain  until  more  extensive  experimental 
data  become  available. 

In  conclusion,  success  achieved  in  this  study  coupled  with  prior 

5 

success  in  predicting  trnasltlon  on  Mach  5  ground-test  and 
Mach  20  flight-test  vehicles^  further  substantiates  the  notion 
that  the  method  can  be  used  as  a  design  tool  for  a  wide  range  of 
aerodynamic/hydrodynamic  vehicles . 
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APPENDIX:  NEAR-SURFACE  BEHAVIOR  OF 

SOLUTIONS  TO  THE  MODEL  EQUATIONS 


This  Appendix  presents  an  analytical  solution  to  the  turbulence 
model  equations  which  Is  valid  very  near  a  solid  boundary.  To 
determine  how  far  from  a  solid  boundary  the  analytical  solution 
applies,  numerical  solutions  to  the  model  equations  have  been 
generated  for  a  variety  of  laminar,  transitional  and  turbulent 
boundary  layers.  As  a  result,  the  required  number  of  mesh  points 
needed  for  accurate  numerical  solution  has  been  reduced  with  an 
attendant  reduction  In  computing  time. 


Computations  with  DCW  Industries'  EDDYBL  computer  code,  which 

embodies  the  turbulence  model  equations,  generally  require  mesh 

points  extremely  close  to  a  solid  boundary  to  insure  accurate 

computation  of  the  turbulent  dissipation  rate,  u.  The  reason  for 

having  very  fine  resolution  near  the  boundary  is  the  singular 

behavior  of  oj  at  very  smooth  surfaces.  For  Incompressible  flow 

20 

over  a  rough  surface,  Saffman  and  Wilcox  have  shown  that 


(1) 


0) 


w 


1  + 


(Al) 


where  w  Is  the  value  at  the  surface.  The  quantity  S  Is  the  unlver 
w 

sal  function  of  surface  roughness  and  mass  Injection  defined  in 
Equations  (19)- (21);  for  = 0  and  for  small  roughness  heights,  k, 

S  is 

s  ‘  (''2) 

Note  that  for  a  perfectly  smooth  surface 


a)-*-y^asy-^0  (A3) 

It  is  therefore  a  y~^  singularity  we  are  forced  to  resolve  In 
ilDDYlH,  computations. 


If  the  analytical  behavior  of  w  is  accurately  represented  by 
Equations  A1  and  A2 ,  we  could  obviate  the  stringent  resolution  prob¬ 
lem  by  using  these  equations  rather  than  computing  co  with  finite 
differences.  Hence,  efforts  focused  first  upon  generalizing  Equa¬ 
tion  A1  for  flows  with  heat  transfer  and  then  testing  the  accuracy 
of  the  formula  for  several  flows. 

Using  singular  perturbation  techniques  we  find  that  the  generalized 
solution  for  flows  with  heat  transfer  Is 

+  —  (A^) 

u 

T 

Where  u  is  velocity.  Note  that  for  Incompressible  flat -plate 
boundary  layers 

u/u^  =  u^y/v  as  u^y/v  -*-0  (A5) 

wherefore  Equation  A4  reduces  to  Equation  Al. 

To  test  the  accuracy  of  Equation  a4,  boundary  layer  computations 
have  been  performed  for  a  variety  of  flows  Including  Incompressible 
and  compressible  FPBL’s,  Incompressible  FPBL  with  uniform  blowing, 
and  flow  over  Body  2  with  heating.  As  shown  In  Figures  Al  and  A2, 
with  the  exception  of  the  blowing  case.  Equation  a4  Is  accurate 
for  values  of  u/u  below  about  4.  For  the  blowing  case,  Equation  A4 

T 

is  valid  to  about  u/u^  =  1. 

Hence,  Equation  A4  can  be  used  for  computing  w  below  u/u^  of  about 
4  for  the  flows  of  Interest  in  this  study  (incompressible,  smooth 
wall).  Rather  than  having  to  place  mesh  points  as  close  to  the  wall 
as  u^y/v =  1/10,  there  Is  no  need  to  have  points  any  closer  than 
u  y/v =  1.  This  results  In  a  substantial  (factor  of  2)  decrease  In 
computing  time. 
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